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Abstract 

We present an exact derivation for the asymptotic large distance 
behavior of the spin two-point correlation function in the XY-model. 
This allows for the exact obtainment of the critical exponent rj — 
1/4 at the Kosterlitz-Thouless transition that occurs in this model 
and in the 2D neutral Coulomb gas and which has been previously 
obtained by scaling arguments. In order to do that, we use the 
language of sine-Gordon theory to obtain a Coulomb Gas description 
of the XY-model spin correlation function, which becomes identified 
with the soliton correlator of that theory. Using a representation 
in terms of bipolar coordinates we obtain an exact expression for 
the asymptotic large distance behavior of the relevant correlator at 
/3'^ = Stt, which corresponds to the Kosterlitz-Thouless transition. 
The result is obtained by approaching this point from the plasma 
(high-temperature) phase of the gas. The vortex correlator of the 
XY-model is also obtained using the same procedure. 

KEY "WORDS: soliton correlation functions, sine-Gordon theory, XY-model, 
2D Coulomb gas, Kosterlitz-Thouless transition 

1 INTRODUCTION 

The sine-Gordon (SG) model is certainly one of the best studied of (1 + 1)- 
dimensional physics. The interest in this field theoretical model has been en- 
hanced by its connections with the two-dimensional (2D) neutral Coulomb gas 
(CG) and also with the 2D XY-magnetic system In this framework, it 
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becomes an useful and powerful tool for the study of a great variety of physical 
properties of these two systems, which in principle admit actual realizations in 
nature. The SG model is integrable in the sense that the spectrum and the S- 
matrix are exactly known Nevertheless, there is a lack of exact results for 
correlation functions, except for a specific value of the coupling parameter 

Many interesting results have been obtained recently, concerning the SG sys- 
tem. The thermodynamic Bethe ansatz has been used for obtaining the free 
energy and specific heat of the system Exact form factors of the soliton op- 
erators and other fields have been derived Density correlation functions 
have been calculated using these form factors Several results concerning the 
thermodynamics of the classical CG have also been obtained Among these, 
we mention the exact free energy for < 4n, in the case of point particles |10| 
and for 47r < < Gvr, in the case of extensive ones Charge and particle 
correlators have been obtained in the low temperature (/3^ > Stt) phase '121. 

In this work we present a computation of the two-point spin correlation func- 
tion of the XY-model using its connection to the SG theory. In particular, we 
exploit the relation existing between the spin operators of the former and the 
soliton creation operators of the latter. Then, we use the CG representation of 
these functions in order to derive an exact series describing the large distance 
behavior of them. This series clearly exhibits two distinct types of asymptotic 
behavior, separated by a critical point at = 8tv. At this point, we obtain 
an exact result for the spin-spin correlator of the XY-model, which exhibits a 
power- law behavior with exponent 77 = i, well-known from the scaling analysis 
of the XY-model This special point, indeed, is the Kosterlitz-Thouless (KT) 
critical point I13| and, in the 2D CG language, it corresponds to the temper- 
ature Tc in which the system undergoes a phase transition from a metallic (or 
plasma) phase, composed of charged particles, into an insulating (or dieletric) 
phase, composed of neutral dipoles (bound pairs of charges) . For < Stt (high- 
temperature), our series exhibits a nontrivial large distance behavior, whereas 
for f]'^ > 8n (low-temperature), its large distance behavior is determined by the 
free theory, as expected. 

The KT critical exponent is usually obtained through scaling arguments, 
which are approximate, in the low-temperature phase 2 14 15, |16| or by heuristic 
arguments in the high-temperature phase |17| . The result presented here is the 
first exact evaluation of the KT critical exponent, performed by approaching the 
critical point from the high-temperature (in CG language) phase. 

In section 2, we review the connection between the SG theory with the 2D 
CG and also with the magnetic XY-model. We pay special attention to the corre- 
spondence between the soliton creation operators in the SG theory and the spin 
operators in the 2D XY-model. This allows us to establish the equivalence be- 
tween their respective correlation functions, which can be written in a convenient 
form with the help of the 2D CG picture. In section 3, we make use of bipolar 
coordinates in order to put the correlation functions in a form that easily allows 
for the obtainment of an exact asymptotic solution at the KT critical point. In 
section 4, we consider the case of the vortex correlation function. Finally, some 
concluding remarks are presented in section 5. Two appendices are included. 
In Appendix A, we review the details of the system of bipolar coordinates em- 



2 



ployed in the calculations, whereas in Appendix B, we demonstrate properties of 
the functions K {\x — y\) and K' {\x — y\) appearing in our final expressions for 
the soliton and the vortex correlation functions, respectively, which are used in 
obtaining the asymptotic behavior of these correlators. 

2 SOLITON AND SPIN OPERATORS 

In this section, we are going to review the connection of the SG theory with the 
2D neutral CG and with the 2D XY-model. We will highlight, in particular, the 
relation of the soliton creation operator of that theory with the spin operator of 
the latter model, as well as the representation of its correlation functions in the 
framework of the classical CG. This will be our starting point for the evaluation 
of correlation functions at the KT point. 

We start from the SG euclidean action, given by 



S — I dxdr 



^d^(t)d^i4> + 2q!o cos/ 



(1) 



It is well-known that the vacuum functional of the theory may be written as 

^ ^ ]™o S W J n '^'^^ ''''P ] 8^ 51 ^'^^ 1" [1^' ~ ^^1' + ^1 f ' (2) 

where Ai = 1 for 1 < i < n and Ai = — 1 for n + 1 < i < 2n and £ is a short- 
distance regulator, which is needed in the case of point particles or, equivalently, 
of a local field theory. The renormalized coupling a is related to the one in Q 
by 

Q = ao(e')^. (3) 
In expression Q we recognize Z as the grand-partition function for the classical 
2D CG of point particles with charges ±1, where a is the fugacity and |^ — -p^^, 
where Tea is the Coulomb gas temperature. In obtaining this relation, we are 
using the convention of 1181 for defining the 2D Coulomb potential. In the case of 
particles with a finite diameter, d, the integration region for the 2n 2:i-integrals 
in ||5J must exclude the regions where \zi — Zj\ < d. 

The theory possesses soliton excitations, bearing a topological charge 



2tt 



dxd^(P. (4) 



At the quantum level, the corresponding states are created by the soliton operator 
of Mandelstam ^ 

fi{x,T) — exp I i ^ / dz(j>(z,T)\. (5) 



As we will see below, this operator has an interesting physical interpretation in 
terms of the associated spin system, namely, the XY-model. 
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The exact soliton mass spectrum is known, as well as the masses of soliton 
bound states (breathers), which occur for 0^ < 47r [H]- Their masses are given by 



M„ =2Msin(^^n) , n = l,2,... <i (6) 

where £ = „ ' and M is the soliton mass. 

The SG/CG system is closely related to a 2D spin system, namely the 0(2) 
ferromagnetic Heisenberg model on a square lattice, or XY-model pi. This is 
described by the hamiltonian 

HxY = - J ^ • , (7) 

where ^11^ 



(8) 



The sum in Q runs over nearest neighbors and T is the temperature (note 
that the temperature in the XY-model, T, is not the same as in the CG, Tcg)- 
Observe that \nif = 1 and 6i is proportional to the angle the spin Hi makes with 
the X-axis. 

When we take the continuum limit of ^ and consider the contributions of 
both spin-waves and vortices to the total energy, we arrive at the hamiltonian [5| 



HxY = 



1 - ^ / 2 J 

-ve -ve + Kcos 2-K\ — 

2 I V r 




(9) 



where we have considered that 9i — > 6{x) in the continuum limit and is a scalar 
field related to 9{x) by 



e{x,y)^ / dzdy^{z,y). (10) 

— OO 

From inUi it follows that 

ve -ve = V4)-V(f) (ii) 

and we see, by comparing ||UJ with Q, that the hamiltonian of the XY-model 
coincides with the euclidean action of the SG model, provided we make the 
identifications: j3 = and k — 2ao. It also follows that the euclidean 

vacuum functional of the SG theory, the grand-partition function of the 2D CG 
and the partition function of the XY-model are all identified. 

In the XY-model, similarly to the case of the CG of extensive particles, we 
have a natural short distance cutoff, namely, the lattice spacing r ,2:. Accord- 
ingly, the integration region for the 2n 2:i-integrals in J^J must exclude, in this 
case, the regions where \zi — Zj \ < r. 

Within this picture, a very interesting and useful connection can be estab- 
lished between the XY-spin thermal correlation functions and the SG quantum 
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soliton correlators. Indeed, in view of 1101 we see that the SG soliton operator 
can be written as 

fi{x, t) = exp \i t) I . (12) 



The relevant XY-correlator, therefore, may be written as 

(n(r-) . n{Q))xY = (M(r-)M^(0))sG, (13) 

where the first expression is the XY-thermal correlator , whereas the second one 
is the SG quantum soliton correlator. 

Within the CG description, the soliton correlator 1131 is given by )15| : 



{M(x)/it(y)) = n 



o 2n 



X exp < — ^ AiAj In l^i - £,1^ + A4arg(£i - y) - arg(2i - f)] > , (14) 

I i7^i = l i = l J 

where is given by (j2j. In the previous expression, we are considering either 
the case of the XY-model or the CG of extensive particles, where, as we have 
mentioned, there is a natural short distance cutoff. Consequently, the integrals in 
Ijl4|l and in Z are defined in the region V{r), in which r < \zi — Zj \ < R, r being 
the lattice spacing r in the former case and the particle diameter d in the latter, 
whereas R is the radius of the system. V{r), therefore defines a finite "volume" 
{V ~ ttR^ since r « R) for the system. We have droped the regulator e from 
the Coulomb potential used in 1141 and in Z, since it is no longer needed in the 
presence of the natural short distance cutoff r. 

Observe that, in 1141 the contribution coming from the soliton operators at 
X and y corresponds, in the CG language, to the interaction of the charges of the 
gas with an external string of electric dipoles orthogonal to it, going from x to 
y, plus the self-interaction of this string, which consists in the first term in 1141 

m 

3 THE SOLITON CORRELATION FUNCTION 

In this section, we are going to obtain a representation of the soliton correlator 
1141 . valid for x ^ y, which will enable us to derive an exact expression for the 
asymptotic large distance behavior of (/i/i^) at the KT point, (3^ = Stt. 

Indeed, for x ^ y we may use the bipolar coordinates described in the Ap- 
pendix A, and rewrite 1141 as 

2. in |._.|4„ 



2^-1 ^2n /"^T n + oc 




— •'o,vir) J-oc,v(.) 4[coshr?, ~ cosC.]2 



-1 2 

X — y\ 



sinh rji sinh rjj 



2 [cosh rji — cos ^i] 2 [cosh rjj — cos ^j] 
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+ 



sin 



sin 



2[coshr7i — cos^i] 2[cosh77j — cos^j] 



2n 
i = l 



(15) 



In this expression, as before, the symbol V{r) expresses the fact that the inte- 
grations must respect the condition that r < \zi — Zj \ < R. In terms of the 
integrals, this implies the following restriction for the expressions between round 
brackets in (I15II . which we call, respectively Uij and l3ij: 



\x — y\>>r 









(16) 



By inspecting dJ, we immediately conclude that the \x — j7|-factors decou- 
ple from the (,i,rii integrals. Elementary combinatorics, taking into account the 
neutrahty of the gas, shows that this factor will appear n(n — 1) times in the nu- 
merator and times in the denominator. Combining with the 4n contribution 
coming from the scale factors of the volume elements, we obtain 



F - y\ 



■E 

n=0 



Cn (1^- 271) |a 



(17) 



In this expression the coefflcients d are given by the summand in 1151 after the 
\x — y|-factors have been removed. In view of the restriction on the integration 
region imposed by l|16^ . we see that these coefficients depend, in general, on 
|a; — y|. However, in the large distance regime, where |a; — y| >> r, we infer from 
the right part of Ijl6|l that the integration region appearing in the expression of 
Cn is restricted by constant bounds and, therefore, do not depend on \x — y\. 
Hence, we conclude that the coefficients C„ are constant in this limit. This fact 
is confirmed from an independent point of view in Appendix B. 

FromEq. lO we see that for 0^ = 8n we have the following exact expression, 
valid for the XY-model and the CG of extensive particles, 

Ki\x-y\) 



(18) 



■ y\ 



As we shall prove in Appendix B,_ftr(|a; — j7|)isa real function with an upper 
bound equal to one, which is saturated in the large distance regime, namely, 

y\) - 



where K ( 



K 



1. This immediately allows us to write 



(M(^)M'(y)) 



^ f ^\\ \S-y\>>r 



1 



(19) 



-y\ 



Eq. 11911 is the well-known result of Kosterlitz and Thouless for the XY-spin 
correlation function at the critical temperature Tkt = t^J (remember the relation 
between the XY-model temperature T and /3, namely, T = ^^5^) Notice 
that for determining the critical exponent of the static and uniform magnetic 
susceptibility, only the large distance behavior of the spin correlator is needed. 
Expression l|19|l provides an exact result for this. 
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We remark, at this point, that only the leading term in the large distance 
behavior of the spin correlator has been obtained when we took the |a; — i/l » r 
limit. In order to get the subleading logarithmic correction obtained previously 
Emg, namely 

, i,-^, \x-y\»r C (In |x — i/l) 5 

{fi{x)fiHy))' ^ \J J^> , (20) 

\x — y\i 

we should consider the next term in the j^r^ expansion. According to the bound 
and asymptotic value obtained for the function K{x) in Appendix B, however, 
we conclude that the constant C in 12011 must be C = (ln2_R)~s such that 

K{\x-y\) ^^''\'^-yf ''-^>^ 1. (21) 

(ln2i?)8 

Finally, observe that the exponent 1/4 comes from the free soliton correlator, 
which is the prefactor in (I15II and 1171 1. This expresses the well-known fact that 
the cos jScj) interaction becomes irrelevant at 0^ = Stt. Equivalently, in XY- 
model language, we would say that the whole contribution to the asymptotic 
behavior of the correlators comes from the spin-wave term. For 0^ > Stt, we still 
have the large distance behavior of the correlators determined by the free theory 
(spin- wave term). This fact can be inferred directly from l|17|l . even though 
some logarithmic corrections could be produced by summing in n. For > 
Stt, a careful analysis of the ultraviolet divergences is required |22l in the local 
case, where any short distance regulators must be removed at some point. This, 
however, has no effect on our result 1171 . where short distance singularities are 
absent due to the presence of natural regulators. 



4 THE VORTEX CORRELATION FUNCTION 

The large distance behavior of the two-point correlation function of the vortex 
creation operators in the XY-model can also be obtained exactly within our 
approach, at the KT point. Here, again, the use of the correspondence with the 
SG theory is quite useful. Indeed, from (|SJ it follows that, in SG language, this 
operator is given by 

(7(a;,r) = exp|i ^ ?i(a;,r)|. (22) 

We can obtain a CG representation for < aa^ > analogous to Ijl5|l . The only 
differences are the prefactor exponent, which is now /J^/Stt and the last term, in 
which is replaced by rji and i by 0^ /A-k. Following the same procedure as for 
the soliton correlation function we obtain, at 0^ = Stt, 

Hx)a^i0!)) = ^^^^ , (23) 
k - y\ 

where the function K' {\x — y\) can be obtained from K (\x — y\) through the 
same replacements described above. In Appendix B we show that K' {\x — y\) 
presents the same bound and asymptotic limit as K (\x ~ y\). 
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Thus we have 

(a(f)at(j7)> l^^-J^^" ^ . (24) 
5 CONCLUDING REMARKS 

An interesting extension of this work, which we are now considering, is the ob- 
tainment of the soliton or XY-vortex four-point correlation functions at the KT 
point, using the same methodology. It would also be interesting to compare the 
large distance behavior of our series (1151 with the exact solution of the associated 
free massive fermion theory at the Luther- Emery point 0^ = 47r |19l I23| . This 
would allow us to determine the coefficients C„ at this point. Another interesting 
issue to be explored would be the large distance behavior of the soliton and vortex 
correlators in the plasma phase ((iP < Svr), where an exponential decay should be 
recovered. Finally, an important case that deserves further investigation is the 
local SG theory or CG of point particles for 0^ > Svr, where a careful treatment 
of the short distance singularities must be performed. 
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APPENDIX A 

Here we give details about the bipolar coordinates, which we use for comput- 
ing the soliton correlation function 11411 . Given the position vector r in the plane 
and two points (poles) at x and y, we define the coordinates r/) as |24| 

^ = arg(r - y) - arg(r - x) 

r? = ln^^, (A.l) 

\r - y\ 

with < ^ < 27r and -co < rj < oo. In terms of these coordinates, the position 
vector is given by 

'"""^l (sinh?7 , sine) (A.2) 



2 [cosh ri — cos ^ 
and the volume element reads 



= .r ^' Tvjd^dr,. (A.3) 

4 [cosh r; - cos e] 2 



APPENDIX B 

From llTHll . llTTl and llTHl we obtain 



n — ^ \ > ) I— I 2—1 



i<7 I x — yl'^ 1 liKj \x — y\ 



K(\x-y 

»l llij- \3-v\^ 

{n n ^ 

-y)- arg(fi - f)] - j^[arg(j7i - y) - arg(j7i - f)] >(B.l) 

In this expression, we went back to usual coordinates and associated Xi and j/i 
with the positive and negative charges, respectively. We also introduced the 
symbols 

[Xi,yj] = \xi~yj\'^ >r'^. (B.2) 

Let us first prove that K is real. This can be easily done by taking the 
complex conjugate of l|B.lfl and subsequently performing the change of variables 
Xi ^ yi. 

Considering now the contribution of the phase factors to the integrals in IB.H 
we can establish the following bound: 



Ki\x~y\)<Z l^T^ {{dx^YYdy. n" \x- v] ^ ^ 



Observe that the [a; — y|-terms in ijB.l^ have completely canceled out. Notice also 
that the expression given by the sum in IIB.3t is nothing but Z and therefore we 
immediately infer that K {\x — y\) < 1, which is the announced bound for the 
function in IjlSjl . 
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We now consider the asymptotic behavior of K, for |a; — j7| >> r. The phases 
in IB. Ill cancel out in the leading order in this limit and we immediately see that 
the above bound is saturated in the leading asymptotic regime, namely 

/, ^ \3-y\>>r 

K{\x~y\) > 1. 

In order to complete our analysis, let us establish now an upper bound for 
Z itself. Using the expression for Z given by the sum in <B.3L we can infer the 
following bound for Z, by making the replacement [xi,yj\ <-> in the n i = j 
terms in the denominator: 



^2n 



z<^ (^) 



^ (n!)2 

71 = 



„ n n 

), (B.4) 



where 



/„(a;i,...,a;„;yi,...,y„) = . (B.5) 

Using the Mean- Value Theorem, we may write the n-th integral in ljB.4^ as 

In = V^'^j:, (B.6) 

where /„ is the average value of /„ and V is the "volume" (~ -kB?) of the system. 
Observe that f\ = 1, whereas /2 is given by 

^ [xi,X2][yi,y2] ,g ^- 

^ [xi,y2][x2,yi] 

and 

h = h{xi,x2\yi,y2) h{xi,xz\yi,yA) h{x2,xz\y2,yi)- (B.8) 
The generalization for higher values of n is straightforward and we have in general 

Consider now l2- Performing the changes of variables x\ y\ and yi xi, 
we readily conclude that 

h=l^^- (B.9) 



This implies that /2 < 1 and also that (/2)^ < 1. 

Since, as we have just seen, /„ is in the form (/2)^ we conclude that /„ < 1 
and, consequently, 7„ < V^". Inserting this bound in IIB.411 . we obtain (notice 
that Z is positive) 



(n!)2 " r2 ) 



< ^ < 2. ^^TTT^ = lo ( — ) , (B.IO) 



where lo{x) is a modified Bessel function of the first kind and we have set V — 
■rB?, since R » r. From this we see that Z is finite for the XY- model and 
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the CG of extensive particles (where there is a natural short distance cutoff r) 
whenever the volume is finite . 

Now let us consider K' {\x — y\). From the explanations given after 12211 . we 
conclude that 



K'{\x~y\)=Z- 



Observe that the last factor in the previous expression tends to one in the 
limit \x — y\ » r. Hence, we immediately conclude, after the \x — i7|-terms are 

canceled out, that K' {\x — y\) ' Zl^^ thus establishing (1241 . 
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